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I. INTRODUCTION 

The study of extrasolar giant planets (EGPs*) and brown dwarfs 
via reflex steUar motion, broad-band photometry, and spectroscopy 
has finally come into its own. Doppler spectroscopy alone has revealed 
about 15 objects in the giant planet/brown dwarf regime, including 
companions to r Boo, 51 Peg, v And, 55 Cnc, p CrB, 70 Vir, 16 Cyg 
B, and 47 UMa (Noyes et al. 1997; Butler et al. 1997; Cochran et al. 
1997; Marcy & Butler 1996; Butler & Marcy 1996; Mayor & Queloz 
1995; Latham et al. 1989). 

The direct detection of G1229 B (Oppenheimer et al. 1995; Naka- 
jima et al. 1995; Matthews et al. 1996; Geballe et al. 1996; Marley 
et al. 1996; AUard et al. 1996; Tsuji et al. 1996) was a milestone since 
G1229 B displays methane spectral features and low surface fluxes that 
are unique to objects with effective temperatures (in this case, Teff~950 
K) below those at the solar-metallicity main sequence edge ('^1750 K, 
Burrows et al. 1993). In addition, the almost complete absence of 
spectral signatures of metal oxides and hydrides (such as TiO, VO, 
FeH, and CaH) is in keeping with theoretical predictions that these 
species are depleted in the atmospheres of all but the youngest (hence, 
hottest) substellar objects and are sequestered in condensed form be- 
low the photosphere (Lunine et al. 1989; Marley et al. 1996). The 
wide range in mass and period, as well as the proximity of many of the 



* We use this shorthand for Extrasolar Giant Planet, but the terms "ex- 
oplanet" or "super-Jupiter" are equally appropriate. 
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planets/brown dwarfs to their primaries, was not anticipated by most 
planetary scientists. Though the technique of Doppler spectroscopy 
used to find these companions selects for massive, nearby objects, their 
variety and existence is a challenge to conventional theory. Since di- 
rect detection is now feasible, and has been demonstrated by the recent 
acquisition of G1229 B, it is crucial for the future of extrasolar giant 
planet searches that the spectra, colors, cvohition, and physical struc- 
ture of objects from Saturn's mass (0.3 Mj ) to 70 Mj be theoretically 
investigated. 

II. EARLY CALCULATIONS OF THE EVOLUTION AND 
STRUCTURE OF EXTRASOLAR GIANT PLANETS 

A. Gray Models 

EGPs radiate in the optical by reflection and in the infrared by 
the thermal emission of both absorbed stellar light and the planet's 
own internal energy. In Burrows et al. (1995) and Saumon et al. 
(1996), the EGPs were assumed to be fully convective at all times. We 
included the c;fTects of "insolation" by a central star of mass and 
considered semi-major axes (a) between 2.5 A.U. and 20 A.U. Giant 
planets may form preferentially near 5 A.U. (Boss 1995), but as the new 
data dramatically afflrm, a broad range of semi-major axes can not be 
excluded. In these calculations, we assumed that the Bond albedo of 
an EGP is that of Jupiter (~0.35). For the Burrows et al. (1995) 
study, we evolved EGPs with masses from 0.3 Mj (the mass of Saturn) 
through 15 Mj . Whether a 15 Mj object is a planet or a brown dwarf 
is largely a semantic issue, though one might distinguish gas giants and 
brown dwarfs by their mode of formation (e.g., in a disk or "directly"). 

If 51 Peg b is a gas giant, its radius is only 1.2 Rj and its luminosity 
is about 3.5 x 10~^ Lq. This bolometric luminosity is more than 1.5 x 
10'' times the present luminosity of Jupiter and only a factor of two 
below that at the edge of the main sequence. The radiative region 
encompasses the outer 0.03% in mass, and 3.5% in radius. The study 
by Guillot et al. (1996) demonstrated that 51 Peg b is well within its 
Roche lobe and is not experiencing significant photoevaporation. Its 
deep potential well ensures that, even so close to its parent, 51 Peg b 
is stable. If 51 Peg b were formed beyond 1 A.U. and moved inward 
on a timescale greater than ~ 10^ years, it would closely follow a iip ~ 
Rj trajectory to its equilibrium position. 

B. Non Gray Models 

However, to credibly estimate the infrared band fluxes and improve 
upon the black body assumption made in Burrows et al. (1995) and 
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Saumon ct al. (1996), wo have recently performed non-gray sim- 
ulations at solar-metallicity of the evolution, spectra, and colors of 
isolated EGP/brown dwarfs down to Tegs of 100 K (Burrows et al. 
1997). Figure 1 portrays the luminosity versus time for objects from 
Saturn's mass (0.3 Mj ) to 0.2 Mq for this model suite. The early 
plateaux between 10^ years and 10^ years are due to deuterium burn- 
ing, where the initial deuterium mass fraction was taken to be 2x10^"''. 
Deuterium burning occurs earlier, is quicker, and is at higher luminos- 
ity for the more massive models, but can take as long as 10^ years for a 
15 Mj object. The mass below which less than 50% of the "primordial" 
deuterium is burnt is ~13 Mj (Burrows et al. 1995). On this figure, 
we have arbitrarily classed as "planets" those objects that do not burn 
deuterium and as "brown dwarfs" those that do burn deuterium, but 
not light hydrogen. While this distinction is physically motivated, we 
do not advocate abandoning the definition based on origin. Neverthe- 
less, the separation into M dwarfs, "brown dwarfs" , and giant "planets" 
is useful for parsing by eye the information in the figure. 

In Fig. 1, the bumps between 10~^ Lq and 10~^ Lq and between 
10* and 10^ years, seen on the cooling curves of objects from 0.03 MqIo 
0.08 Mq, are due to silicate and iron grain formation. These effects, 
first pointed out by Lunine et al. (1989), occur for Tefjs between 2500 
K and 1300 K. The presence of grains affects the precise mass and 
luminosity at the edge of the main sequence. Since grain and cloud 
models are problematic, there still remains much to learn concerning 
their role and how to model them (Lunine et al. 1989; AUard et al. 
1997; Tsuji et al. 1996). 

III. NEW INSIGHTS 

The studies of Burrows et al. (1997) and Marley et al. (1996) 

reveal major new aspects of EGP/brown dwarf atmospheres that bear 
listing and that uniquely characterize them. Below Tegs of 1300 K, the 
dominant equilibrium carbon molecule is CH4, not CO, and below 600 
K the dominant nitrogen molecule is NH3, not N2 (Fegley & Lodders 
1996). In objects with Teg < 1300 K, the major opacity sources are H2, 
H2O, CH4, and NH3. For Tefjs below ~400 K, water clouds form at or 
above the photosphere (defined where T = Tcff ) and for Ton s below 200 
K, ammonia clouds form (viz., Jupiter). Collision-induced absorption 
of H2 partially suppresses emissions longward of ~10 /um. The holes in 
the opacity spectrum of H2O that define the classic telluric IR bands 
also regulate much of the emission from EGP/brown dwarfs in the near 
infrared. Importantly, the windows in H2O and the suppression by H2 
conspire to force flux to the blue for a given Tcff. The upshot is an 
exotic spectrum enhanced relative to the black body value in the J 
and H bands (-^1.2 /xm and ~1.6 /xm, respectively) by as much as two 
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to ten orders of magnitude, depending upon Tdf. The enhancement 
at 5 /xm for a 1 Gyr old, 1 Mj extrasolar planet is by four orders of 
magnitude. As Teg decreases below ~1000 K, the flirx in the M band 
(~5 fim) is progressively enhanced relative to the black body value. 
While at 1000 K there is no enhancement, at 200 K it is near 10^. 
The J, H, and M bands are the premier bands in which to search for 
cold substcUar objects. The Z band (^^1.05 /im) is also in excess of 
the black body value over this T^g range. Even though K band (~2.2 
/xm) fluxes are generally higher than black body values, H2 and CH4 
absorption features in the K band decrease its importance relative to 
J and H. As a consequence of the increase of atmospheric pressure 
with decreasing Teff , the anomalously blue J — K and H — K colors get 
bluer, not redder. The K and J versus J — K infrared H R diagrams 
loop back to the blue below the edge of the main sequence and are 
not continuations of the M dwarf sequence into the red. The difference 
between the black body curves and the model curves is between 3 and 
10 magnitudes for J versus J — K, more for K versus J — K. G1229 B 
flts nicely among these theoretical isochrones. The suppression of K by 
H2 and CH4 features is largely responsible for this anomalous blueward 
trend with decreasing mass and Teg. 

IV. ATMOSPHERIC COMPOSITIONS OF EGPS, CON- 
DENSATION, AND CLOUDS 

The molecular compositions of the exotic, low-ionization atmo- 
spheres of EGPs and brown dwarfs can serve as diagnostics of temper- 
ature, mass, and elemental abundance and can help define a spectral 
sequence, just as the presence or absence of spectral features associated 
with various ionization states of dominant, or spectroscopically active, 
atoms and simple molecules does for M through O stars. However, the 
multiplicity of molecules that appear in their atmospheres lends an ad- 
ditional complexity to the study of substellar mass objects that is both 
helpfully diagnostic and confusing. Nowhere is the latter more appar- 
ent than in the appearance at low temperatures of refractory grains 
and clouds. These condensed species can contribute significant opac- 
ity and can alter an atmosphere's temperature/pressure profile and its 
albedo. Grain and cloud droplet opacities depend upon the particle size 
and shape distribution and these are intertwined with the meteorology 
(convection) in complex ways. Furthermore, condensed species can 
rain out and deplete the upper atmosphere of heavy elements, thereby 
changing the composition and the observed spectrum. In brown dwarf 
and EGP atmospheres, abundance and temperature/pressure profiles, 
particle properties, spectra, and meteorology are inextricably linked. 

The formation of refractory silicate grains below 2500 K was al- 
ready shown by Lunine et al. (1989) and Burrows et al. (1989) to in- 
fluence the evolution of late M dwarfs and young brown dwarfs through 
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their "Mic" opacity. The blanketing effect they provide lowers the effec- 
tive temperature (Teg) and luminosity (L) of the main sequence edge 
mass from about 2000 K and 10"^ Lq to about 1750 K and 6 x 10"^ 
Lq , an effect recently verified by Chabricr et al. (1998). In addi- 
tion, grain opacity slightly delays the cooling of older brown dwarfs, 
imprinting a slight bump on their luminosity /age trajectories (see Fig. 
1). The presence of grains in late M dwarf spectra was invoked to ex- 
plain the weakening of the TiO bands and the shallowing of their H2O 
troughs in the near infrared (Tsuji et al. 1996; Jones & Tsuji 1997). 
Tsuji and collaborators concluded that titanium was being depleted 
into refractories, a conclusion with which we agree. 

With a Teff of ~950 K, a luminosity below 10~^ Lq , and spectra 
or photometry from the R band through 5 /im, G1229 B hints at or 
exemplifies all of the unique characteristics of the family: metal (Fe, 
Ti, V, Ca, Mg, Al, Si) depletions, the dominance of H2O vapor, the 
appearance of CH4 and alkali metals, and the signatures of clouds. 
Clouds of low-temperature condensible species above the photosphere 
are the most natural explanation for the steep drop below 1 /xm in 
the Keck spectra between 0.83 /im and 1 jim (Oppenheimer et al. 
1998). These clouds might not be made up of the classic silicate re- 
fractories formed at much higher temperatures, since these species may 
have rained out. From simple Mic theory, their mean particle size must 
be small (~ 0.1 /um) in order to influence the "optical" without much 
perturbing the near infrared. In addition, such a population of small 
droplets can help explain why G1229 B's near-infrared troughs at 1.8 
/um and 3.0 /xm are not as deep as theory would otherwise have pre- 
dicted. Just as Tsuji and collaborators have shown that silicate grains 
at higher temperatures can shallow out the II2O troughs, so to can 
species that condense at lower temperatures (< 1000 K ?) explain the 
shallower-than-predicted G1229 B H2O troughs. What those species 
might be can be illuminated by chemical abundance studies (Burrows & 
Sharp 1998). Note that a cloud grammage in these small-radius, low- 
temperature refractories of only ~ 10~^ gm cm~^ would be adequate 
to explain the anomalies. 

A. Abundance and Condensate Results 

At temperatures of ^^1500 2000 K, the standard refractories, such 
as the silicates, spinel, and iron, condense out into grain clouds which 
by their large opacity lower the Tes and luminosity of the main se- 
quence edge (Burrows et al. 1993; Burrows & Sharp 1998) and alter 
in detectable ways the spectra of objects around the transition mass 
(Jones & Tsuji 1997). As T^g decreases below that at the stellar edge, 
the classical refractories are buried progressively deeper below the pho- 
tosphere and less refractory condensates and gas-phase molecules come 
to dominate (Marley et al. 1996; Burrows et al. 1997). 
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Below temperatures of ^1500 K, the calculations of Burrows & 
Sharp (1998) demonstrate and confirm that the alkali metals, which 
are not as refractory as Fe, Al, Ca, Ti, V, and Mg, emerge as im- 
portant atmospheric and spectral constituents. At still lower temper- 
atures, chlorides and sulfides appear, some of which will condense in 
the cooler upper atmosphere and form clouds that will affect emergent 
spectra and albedos. Cloud decks of many different compositions at 
many different temperature levels are expected, depending upon Tefj 
(and weakly upon gravity). Clouds of chlorides and sulfides at temper- 
ature levels below ^1000 K may be responsible for the steeper slope 
observed in the spectrum of G1229 B at the shorter wavelengths (Op- 
penheimer et al. 1998). At slightly higher temperatures, MnS, ZnS, 
NaAlSisOs, KAlSisOs, V2O3, and MgTi205 may play a role, but only 
if their constituents are not scavenged into more refractory compounds 
and rained out deeper down. 

As Toff decreases (either as a given mass cools or, for a given age, 
as we study objects with lower masses), the major atmospheric con- 
stituents of brown dwarfs and EGPs change. This change is reflected 
in which spectral features are most prominent and in the albedos of sub- 
stellar objects near their primaries. Equilibrium chemical sequences are 
predominantly a function of temperature and can help to define a spec- 
tral sequence for substellar objects from the main sequence edge near 
2000 K to EGPs with Tegs of a few hundred Kelvin. The appearance 
and disappearance of various molecules and refractories delineates an 
effective temperature sequence and the new proposed "L" dwarf spec- 
tral classification (Kirkpatrick et al. 1998) may correspond to a subset 
of such a compositional sequence. Very crudely, the "L" spectral type 
suggested by Kirkpatrick et al. (1998) would correspond to Tons be- 
tween about ~1500 K and ~2200 K. All but the youngest and most 
massive brown dwarfs and only the very youngest EGPs could have 
this proposed spectral designation. Most brown dwarfs and EGPs will 
be of an even later spectral type, yet to be coined, a spectral type that 
would include G1229 B. 

B. Clouds 

Cloud formation refers to the production of solid or liquid particles 

in the atmosphere of a brown dwarf or extrasolar planet. Cloud forma- 
tion in the Earth's troposphere is a complex phenomenon involving a 
single species (water), two phases, multiple particle size distributions, 
incompletely-known particle properties and myriad cloud morpholo- 
gies. The range of radiative properties of terrestrial clouds is com- 
mensurately broad, and different types of clouds can serve to either 
warm or cool the surface of the Earth. Much of the complexity lies 
in the heterogeneous distribution of such clouds in vertical and hori- 
zontal directions; the problem of treatment of the radiative transfer of 
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broken clouds remains a difficult one (Goody & Yung 1989). Further- 
more, cloud formation itself contributes latent heat to an atmosphere 
and, hence, can destabilize a radiative region, requiring treatment of 
the heat transfer by moist convection (Emanuel 1994). Cloud-forming 
species are known or suspected to be present in eight of the nine plan- 
ets of the solar system, as well as two of the moons of the outer solar 
system (Titan and Triton). 

The giant planets of our solar system remain the best guide to the 
study of clouds in brown dwarfs and extrasolar giant planets, because 
the background gas is hydrogen helium in all cases, and the atmo- 
spheres are not thin layers atop a solid surface. Principal cloud forming 
species in the upper Jovian and Saturnian atmosphere are ammonia, 
water, and possibly ammonium hydrosulfide: in Uranus and Neptune 
they are methane and either ammonium hydrosulfide or hydrogen sul- 
fide (Baines et al. 1995). However, a number of minor species may 
participate in cloud formation, and stratospheric photochemical hazes 
are generated by the action of solar ultraviolet radiation. In situ obser- 
vations of Jupiter's cloud structure by the Galileo entry probe revealed 
the extent of heterogeneity across the surface of a gas giant planet; the 
probe failed to detect anything but very a tenuous and narrow cloud 
layer (Ragent et al. 1995). Since ample visual and spectroscopic evi- 
dence exist for widespread clouds on Jupiter (Carlson et al. 1995), the 
Galileo probe must have found a region of unusual dryness. This notion 
is supported by the observation that the probe fell into a "hot spot," 
a region of imusually high 5 -micron emission (Orton et al. 1996). 

Modeling a potentially broad range of compositions and properties 
from the small number of detailed observations of clouds available to 
date is daunting. One strategy is to consider end-member cloud mod- 
els, characterized by simple condensation in a radiative atmosphere 
on the one hand, and by large scale vertical motions driven by ther- 
mal convection or other processes on the other (Lunine et al. 1986; 
1989). Such models can be reduced to a few parameters which can then 
be constrained, albeit weakly, with spectral observations of the atmo- 
sphere. In the radiative case the cloud mass density versus altitude is 
determined by the saturation vapor pressure and an assumed value for 
the meteorological supersaturation, which determines the condensable 
available to form clouds, but is determined for the Earth empirically 
(Emanuel 1994). This factor is difficult to determine a priori for an ob- 
ject in which the clouds cannot be directly studied, but can be crudely 
estimated (Rossow 1978). A simple prescription for the cloud number 
density for a convectively unstable brown dwarf atmosphere is given 
by Lunine et al. (1989). In each case the mean particle size can be 
determined semi- analytically (Rossow 1978) or by use of a numerical 
particle growth model (Yair et al. 1995). The problem with the former 
approach is that the mean particle size can be significantly underesti- 
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mated, based on terreKtrial experienee, and a size distribution is not 
obtained. The numerical approach suffers from the problem of poor 
specificity of input parameters when applied to extrasolar planets and 
brown dwarfs. 

Spectroscopy of Gliese 229 B suggests that grain formation may be 
occurring in its atmosphere. Grains both alter the spectral contrast of 
molecular lines and, through condensation, remove gas phase species 
that may be directly responsible for various absorption features. Be- 
cause of condensation, silicate or iron features in the spectrum may 
disappear at around or above 1000 K, with gaseous water bands doing 
the same around 400 K. The effects of clouds on the albedo or reflec- 
tivity of a brown dwarf or extrasolar planet are also large. Modeling 
by Marley et al. (1998) suggests large variations (factors of two or 
more) in the albedo of a brown dwarf or extrasolar planet, depending 
upon the type of cloud (large-particle convective versus small-particle 
laminar clouds). The four giant planets all have surface reflectivities in- 
fluenced by a combination of cloud and gas opacity sources, to varying 
extents. 

Further progress in quantifying the extent of the effects of clouds 
on brown dwarf atmospheres will require more realistic models of cloud 
formation, transport and growth of particles, data on indices of re- 
fraction of cloud particles, and incorporation of moist convective insta- 
bility into radiative-convective models. These models will require ad- 
ditional computational resources beyond the already substantial speed 
and memory requirements of fully non gray radiative convective codes. 
Perhaps most challenging from a mathematical point of view is to char- 
acterize and incorporate the size-frequency spectrum of broken clouds. 
The problem is somewhat akin to the classical one of computing Rosse- 
land mean absorption coefficients in gray atmosphere models, where 
the greatest flux contribution comes from the frequency windows of 
smallest opacity. Qiiantifying a broken cloud atmosphere through a 
single number characterizing percentage of area covered by clouds is 
likely to be highly inaccurate, because the amount of thermal radiation 
escaping will be highly sensitive to the gaps between the clouds. In 
the near-term, the most promising observations for constraining such 
a complex situation will be high spectral resolution data obtaining from 
large telescopes with sensitive electronic detectors. However, the the- 
oretical challenge of accurately modeling broken clouds with existing 
computational resources remains unsolved. 

V. MODELS OF THE BROWN DWARF GLIESE 229 B 

The wide gap between stars and brown dwarfs near the edge of the 
main sequence on the one hand, and Jupiter on the other hand, is for- 
tuitously occupied by the cool brown dwarf, Gl 229 B. It is suflaciently 



9 



bright to allow high resolution and high signal-to-noisc spectroscopy. 
This fascinating object displays phenomena common to giant planets 
and to very-low-mass stars and represents a benchmark for modeling 
atmospheres of EGPs. Model spectra for G1229 B (Marlcy ct al. 1996; 
AUard et al. 1996; Tsuji et al. 1996) reproduce the overall energy dis- 
tribution fairly well and all agree that 1) Tes ~ 950 K, 2) the silicate 
opacity is small compared to the gaseous molecular opacity and can 
be ignored in a first approximation, and 3) the gravity of G1229 B is 
poorly constrained at present. The models, however, fail to reproduce 
the visible flux, the observed depth of the strongest molecular absorp- 
tion bands, as well as the detailed structure of the observed spectrum. 
The calculation of better models for G1229 B is currently limited by 
the inadequate knowledge of CH4 opacities and a very limited under- 
standing of the role of grain opacity in such a cool object. 

Aside from possible grain absorption and scattering, the spectrum 
of G1229 B is shaped entirely by absorption by H2, H2O, CH4 and, 
in the mid-infrared, NH3. While the opacities of H2 and H2O are 
now well understood over the entire temperature range of interest for 
G1229 B (500 T 2000 K), the current databases for CH4 and 
NH3 are limited to T < 300 K (HITRAN database, Rothman et al. 

1992) and have an incomplete wavelength coverage. In the case of 
CH4, this incompleteness is partly remedied by complementing the line 
list with frequency-averaged opacity (Karkoschka 1994; Strong et al. 

1993) . Because CH4 plays an important role in shaping the 1 to 4/zm 
spectrum of G1229 B, the near- infrared fluxes of the present models are 
of limited accuracy, precisely where G1229 B is brighter and most easily 
observable {Z, J, H, and K bands). 

The models of Allard et al. (1996) and Tsuji et al. (1996) include 
TiO and VO opacities (as well as metal hydrides) which have very 
strong bands in the red part of the spectrum. As a consequence, their 
synthetic spectra reproduce the rapid decrease of the flux shortward of 
1/um, as evidenced by R and / photometry (Matthews et al. 1996), 
quite nicely. However, the visible spectra of Schultz et al. (1998) and 
Oppcnheimcr et al. (1998) show a strong H2O band but not the bands 
of oxides and hydrides predicted by those models. As anticipated by 
Marley et al. (1996), refractory elements condense in the atmosphere 
of G1229 B and arc removed from the gas phase (Burrows & Sharp 
1998) and a significant fraction of the condensed particles may settle 
to unobservable depths in the atmosphere. Fluxes from models which 
do not include grain opacity in the visible predict visible fluxes which 
are grossly in error. 

The lack of any other important molecular spectral features in the 
visible (Griffith et al. 1998; Saumon et al. 1998) suggests that a source 
of continuum opacity is required. An opacity source that fits these 
requirements is sub- micron grains. Mie scattering theory predicts that 
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0.1 micron radius grains can provide substantial opacity below about 1 
micron yet still be transparent at longer wavelengths (Fig. 2). In fact, 
such behavior is commonly observed in the solar system and has been 
the solution for such diverse puzzles as the radar scattering behavior 
of Saturn's rings particles and the spectrum of Venus' atmosphere. 

Figure 2 demonstrates the dependence of the opacity of small parti- 
cles upon their size. Submicron particles interact strongly with optical 
radiation, while affecting only slightly longer wavelengths. Griffith et 
al. (1998) investigated a range of possible particle optical properties 
and sizes and found that very red, abosorbing, submicron particles can 
indeed lower the optical flux of Gliese 229B, while only slightly affecting 
the near-infrared spectrum. Griffith et al. speculate that since G1229 
B receives twice the UV flux of Titan, an object whose atmosphere 
is dominated by a photochemical smog of condensates, the grains in 
the brown dwarf's atmosphere could also consist of photochemically- 
derivcd non- equilibrium species. Such carbonaceous material is seen 
throughout the solar system and is generally recognizeable by its very 
red color. 

A. Gravity 

It is highly desirable to constrain the value of the surface gravity 
of G1229 B to an astrophysically useful range. As reported by Allard 
et al. (1996), the spectral energy distribution of G1229 B models is fairly 
sensitive to the gravity. The most gravity-sensitive colors are H — K 
and J-i7with A(F-ii')/A(log.g) = -0.39 and A(J-ii:)/A(log.g) = 
—0.26, respectively. However, the uncertainties in the photometry of 
Matthews et al. (1996) in these color indices are comparable to this 
gravity sensitivity. The present uncertainties in CH4 opacities and in 
the role of grains also limit the ability of the models to predict reliable 
near-infrared colors. 

Nevertheless, it is possible to constrain the gravity by analyzing 
the spectrum in the 1.9 - 2.1 /xm region (Saumon et al. 1998). In 
this unique part of the spectrum of G1229 B, the absorption is com- 
pletely dominated by two well- understood opacity sources: H2O and 
the collision-induced absorption by H2. Methane absorption is more 
than two orders of magnitude weaker and unlikely to become signifi- 
cant even when high-temperature CH4 opacities become available. The 
importance of H2O opacity in this region is confirmed by the remark- 
able correspondance of the features of the observed spectrum and of 
the opacity of water (Geballe et al. 1996). The detailed features of the 
synthetic spectrum are therefore far more reliable in the 1.9 - 2.1 /xm 
region than in any other part of the spectrum. Figure 3 shows how the 
gravity affects the modeled features of water in this narrow spectral 
range. The models shown all have the bolometric luminosity of G1229 B 
given by Matthews et al. 1996. The top curve shows the Teg = 870 K, 
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log g = 4.5 model. To emphasize the effects of changing the gravity, 
the flux ratio between this model and the T^s = 940 K, log 5 = 5 model 
(middle curve) and the Teg = 1030 K, logs' = 5.5) model (lower curve) 
are shown. If there were no gravity sensitivity, the two lower curves 
would be flat. This region contains thousands of H2O features and 
their number and depth increases with the spectral resolution. Corre- 
spondingly, the gravity sensitivity washes out at low resolution. High 
signal-to-noise spectroscopy at a resolution of i? ^ 1500 should show 
this effect very well and has the potential of reducing the uncertainty 
in the surface gravity of G1229 B to ±0.25 dex. 

B. Metallicity 

The metallicity of late M dwarfs is notoriously difficult to deter- 
mine and modern efforts are still limited to a small number of stars 
(see, for example, Leggett et al. 1996, Schweitzer et al. 1996, Viti 
et al. 1997). The metallicity of the primary star, G1229 A, appears 
to be approximately solar, but is rather uncertain. Published values 
give [Fe/H] = -0.2 (Schiavon, Barbuy, & Singh 1997) and [M/H]=0.2 
(Mould 1978). The case of G1229 B is even more problematic. The 
relative success of the current models indicate that G1229 B is also ap- 
proximately solar in composition. However, if the brown dwarf formed 
from a dissipativc accretion disk, it may be analogous to the giant plan- 
ets of the solar system and be enriched in heavy elements. On the other 
hand, phase separations in the interior may deplete the atmosphere of 
its heaviest elements. 

Grain formation complicates the definition of the metallicity of 
an atmosphere. Elements are selectively removed from the gas phase 
and introduced in condensed phases whose composition is far more 
difficult to establish by spectroscopy than that of the gas. Since the 
gas-phase abundance of refractory elements is quite sensitive to the 
physics of condensation, it is best to focus on the abundant metals C, 
N, and O which are not significantly depleted by condensed phases in 
the atmosphere of G1229 B. 

The limitations of the CH4 opacity database currently prevents a 
reliable determination of [C/H] in G1229 B. While the opacity of H2O, 
the main oxygen-bearing molecule, is now well understood, the metal- 
licity dependence of the synthetic spectra is muddled by the presence of 
dust affecting the near infrared. The need to untangle the veiling due to 
dust and the effects of metallicity on the H2O absorption bands make it 
very difficult to determine the [0/H] ratio in G1229 B, although Griffith 
et al. (1998) find that a best fit to the optical spectrum of G1229 B is 
achieved with a subsolar oxygen abundance ([0/H] = —0.2). 

Our model spectra predict a strong feature of NH3 near 10.5 fim 
which has not yet been detected by observers (Marley et al. 1996, 
Saumon et al. 1998). The identification of NH3 in the spectrum of 
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G1229 B woiild represent the first detection of this molecule in a com- 
pact object outside of the solar system. It also offers a good possibility 
of measuring the [N/H] ratio, since dust opacity is negligible in G1229 
B at this long wavelength. 

The model atmospheres of Burrows et al. (1997) all assumed a 
solar abundance of the elements. In fact, the metallicity of the sun is 
somewhat higher than that of the average star and it is appropriate 
to consider a greater variety of atmospheric metallicities. We have 
begun this process by computing atmosphere models for brown dwarfs 
in the mass and temperature range of Gliese 229 B. These exploratory 
models were constructed by varying the mixing ratio of all molecules 
uniformly away from that predicted for a solar mixture of the elements 
in thermochemical equilibrium. In actuality, the relative mixing ratios 
of the molecules will not change uniformly as the overall metallicity 
changes. However, such departures will be slight compared to the large 
range in overall metallicity that we have considered. 

Low-resolution spectra for three models, each with T^g = 1000 K 
and g = lOOOmsec"^, but with varying metallicities, are presented in 
Fig. 4. The models were computed following the procedures in Burrows 
et al. (1997), although the treatment of Rayleigh scattering has been 
improved. In the figure, the metallicities are varied over the exception- 
ally large range of [Fc/H] = —2 to 1 to demonstrate the overall trends. 
Generally, as the metallicity decreases, the temperature profile adjusts. 
More flux emerges in the 3 to 5 /zm band as the continuum molecular 
opacity falls. Surprisingly, this redistribution of flux results in a de- 
crease in the flux emerging from the depths of the near-infrared water 
bands. Instead, the lower molecular opacity in the window regions al- 
lows more flux to emerge from deeper, hotter layers of the atmosphere 
resulting in a larger flux in the windows, a cooler upper atmosphere, and 
less flux in the depths of the near-infrared water and methane bands. 
Beyond about 2 /im, the overall flux rises as the metallicity increases. 
The larger metallicity closes off the near-infrared windows, raising the 
upper atmosphere temperature, and increasing the continuum flux. 

The changes in metallicity considered above produce relatively few 
changes in the color differences of G1229 B-like models, particularly 
for solar and subsolar abundances. Significant changes are found only 
for increases in metallicity above solar. J — H and J — K arc most 
sensitive to such metallicity variations, as both colors become redder 
as the metallicity increases above solar. H—K and K—L' are relatively 
insensitive. 

C. Convection 

The molecules found in the atmosphere of a brown dwarf or EGP 
constrain atmospheric structure, dynamics, and chemistry. By identi- 
fying the atmospheric composition, spectroscopy provides information 
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on the physical processes which govern the atmosphere. 

Departures of atmospheric composition from equilibrium are espe- 
cially interesting. CO, PH3, GeH4, and AsHs have all been detected in 
Jupiter's atmosphere at abundances many orders of magnitude higher 
than expected from equilibrium chemistry (see review by Fegley & Lod- 
ders 1994). The presence of these non-equilibrium molecules is taken 
to be evidence of convection. Since convective timescales are shorter 
than chemical equilibrium timescales, these molecules can be dredged 
up from deeper in Jupiter's interior and transported to the visible at- 
mosphere. 

In G1229 B, the detection of CO in abundances in excess of that 
predicted for chemical equilibrium (Noll et al. 1997) implies that the 
visible atmosphere (near 800 to 1400 K) must also be convective. Yet 
many atmosphere models find that the radiative-convective boundary 
lies far deeper, below 1700 K. However, the models of Marley et al. 
(1996) and Burrows et al. (1997) predict an additional, detached, upper 
convection zone. Such a zone would transport CO to the visible atmo- 
sphere from depths where it is more abundant. Other molecules may 
also trace convection, including PH3. The chemical equilibrium profile 
of Cs is very similar to that of CO (Burrows & Sharp 1998). Thus, the 
same convection that dredges CO must also be dredging Cs. However, 
the lack of TiO and other refractory diatomics in the spectrum suggest 
that the atmosphere is not fully convective to the depth (below 2000 
K) where these molecules condense. Taken together, these results may 
support the presence of a detached convection zone. Thus, CO and Cs 
may be tracing the vertical convective structure of the brown dwarf. 
A similar radiative zone, lying below Jupiter's visible turbulent atmo- 
sphere, has been predicted by Guillot et al. (1994). A confirmation of 
such a zone at G1229 B would strengthen the argument for such a zone 
in Jupiter. PH3 is also potentially detectable in G1229 B by space- 
based platforms and will also act as a tracer of convection (Noll et al. 
1997). Measurements of the abundances of this suite of molecules in a 
variety of objects may map out the atmospheric dynamics of substellar 
objects. 

Photochemistry driven by incident radiation can also produce im- 
portant non-equilibrium species. Thus, many hydrocarbons are found 
in the atmospheres of the solar Jovian planets, including C2H2 and 
C2H6 that would not otherwise be expected. A rich variety of photo- 
chemical products will likely be found in the atmospheres of the ex- 
trasolar planets, particularly those with warm atmospheres and large 
incident fluxes. Hazes produced by the condensation of some species 
can produce signatures in the spectra of these objects far in excess of 
what might be expected given their small mixing ratios. 
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VI. ALBEDOS AND THE REFLECTIVITY OF EGPS 

Both scattered light from the primary star and thermal emission by 
a planet contribute to a planet's spectrum. These two components can 
be crudely modeled as the sum of the reflection of a high temperature 
Planck function characterizing light from the primary plus a second, 
lower temperature, Planck function representing the thermal emission 
from the planet itself (Saumon et al. 1996). For the planets of our solar 
system, the two Planck functions are well separated in wavelength. This 
can be seen from the Wien displacement law, AT = 0.29cmdeg. For a 
6000 K primary and a planet radiating at 200 K, the Planck functions 
of the primary and planet peak at 0.48 /zm and 14.5 /zm respectively. 
This separation in the bulk of the radiation field from the planets and 
the Sun has led to a specialized nomenclature in which "solar" and 
"planetary" radiation are often treated separately. 

However, for an arbitrary planet orbiting at an arbitrary distance 
from its primary, there can be substantial overlap of the two Planck 
functions. A general theory of extrasolar atmospheres must consistently 
compute the absorbed and scattered incident light. Marley (1998) has 
generated exploratory EGP atmosphere models that include deposition 
of incident radiation. He finds that in typical EGP atmospheres absorp- 
tion in the strong near-IR water and methane bands produces tempera- 
ture inversions above the tropopausc, similar to Jupiter's stratosphere. 
Generation of a comprehensive suite of model reflected and emitted 
spectra will require that a large range of primary stellar types, orbital 
distances, and planetary masses and ages be investigated. 

Planets are brightest near the peak in the solar Planck function 
and the reflected flux falls off at shorter and longer wavelengths. To 
remove the effects of the solar spectrum and more clearly understand 
the processes acting in the planet's atmosphere, the reflected spectra 
of planets are commonly presented as geometric albedo spectra. The 
geometric albedo is essentially the planetary spectrum divided by the 
solar spectrum. Formally, it is the ratio of the flux received at Earth 
at opposition to the flux that would be received by a Lambert disk of 
the same sizc^ as the planet at its distance from the sun. 

Reflected planetary spectra can also be approximated by com- 
puting wavelength-dependent geometric albedos of EGP atmospheres. 
The resultant planetary spectrum is then the sum of the emitted flux 
plus the product of the incident radiation times the geometric albedo 
with a phase correction. Figure 5 compares the geometric albedo spec- 
tra of Jupiter and Uranus. A purely Rayleigh scattering planet would 
have a geometric albedo of 0.75 at all wavelengths. 

As Fig. 5 demonstrates, planets are not gray reflectors. They 
reflect best near 0.5 fim where Rayleigh scattering dominates the re- 
flected flux. At shorter wavelengths, Raman scattering, which shifts 
some UV photons to longer wavelengths, and absorption by the ubiqui- 
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tous high ahitudc haze found throughout the outer solar system, lower 
the geometric albedo. Longward of 0.6 fim the strength of molecular 
rotational- vibrational bands increases and molecular absorption, rather 
than scattering, begins to dominate the spectrum. Methane and the 
pressure-induced bands of hydrogen are the most important absorbers. 
In between the molecular bands, solar photons reach bright cloud decks 
and are scattered. So the planets remain bright in some band passes. 

Marley et al. (1998) computed geometric albedo spectra for a large 
variety of EGPs, ranging from planets of less than 1 Jupiter mass to 
the most massive brown dwarfs. They considered objects with effective 
temperatures between 100 and 1200 K and found that the UV and op- 
tical spectra of extrasolar giant planets are generally similar to those of 
the solar giant planets. At longer wavelengths, however, the reflected 
flux depends critically on the presence or absence of atmospheric con- 
densates. When condensates are present, photons have the opportunity 
to scatter before they are absorbed. 

The most important condensate in EGP atmospheres is water. Wa- 
ter clouds appear as EGPs cool through effective temperatures of about 
400 K. The sudden appearance of water clouds brightens the planets in 
reflected red and IR light, as shown in Figure 6. This figure presents 
a computed spectrum for an extrasolar planet orbiting at 1 A.U. from 
a G2V primary. Also shown is the emitted flux for an object with 
Teff = 400 K. It is apparent from this plot that the thermal emis- 
sion dominates in the infrared for cloud-free objects. However, once 
clouds form they will both attenuate the emitted flux and reflect a far 
larger proportion of the incident radiation. Of course, there will be a 
continuum such that for objects with the same effective temperature, 
but which are closer to their primary, the reflected flux will surpass 
the emitted flux. The emitted flux will continue to be important at 
lower temperatures for objects further from their primaries. Again, 
individual models arc required for each specific case. 

The ratio of the total reflected light to the total light incident upon 
a planetary atmosphere (Bond albedo) depends sensitively on the spec- 
tral distribution of the incident radiation. For example, most of the flux 
of an A star emerges in the UV and blue. When incident upon a giant 
planet, most of this light may be Rayleigh scattered, resulting in a large 
Bond albedo. In contrast, the predominantly red and infrared photons 
from an M star are far more likely to be absorbed. Thus, the Bond 
albedo of the same planet, when illuminated by two different stars, 
could vary by up to an order of magnitude, as the stellar type of the 
primary varies (Marley et al. 1998). Note that nonequilibrium photo- 
chemical hazes can darken the planets in the UV, further complicating 
the reflected spectra and energy budgets. 
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VII. THE RELATIONSHIP BETWEEN GIANT PLANETS 
OF THE SOLAR SYSTEM AND EGPS 

The advent of the science of EGPs and brown dwarfs, with its many 

new discoveries outside the solar system, should not cause us to lose 
sight of the central role of our local giant planets as exemplars that can 
be studied in unrivaled detail. There are still many questions that re- 
main concerning Jupiter's and Saturn's formation and cooling histories 
that we must answer before we can confidently tackle their analogs be- 
yond the solar system. One of the issues that still surrounds the study 
of our gas giants is the character, radial distribution, and origin of their 
compositions. Do they have a common formation mechanism with the 
newly discovered substellar objects, or is there a distinct mass-related 
boundary between the formation of such giant planets O.OOIMq) 
and brown dwarfs (~ O.OIMq) (see M. Mayor et al., this volume) that 
is reflected in their compositions? Jupiter and Saturn provide unique 
laboratories by which we can address these and related questions of 
fundamental interest to planetary science and astronomy. 

A. Non Solar Metallicities and Implications for Formation 

The solar system exhibits a progressive change in the bulk com- 
position of its four giant planets as a function of heliocentric dis- 
tance. According to interior models constrained by improvements in 
the hydrogen-helium equation of state and by the Galileo entry probe 
results for the abundance of helium in the Jovian atmosphere, Jupiter is 
close to bulk solar composition, with a probable enhancement of C, N, 
and O by about a factor of ~ 5, distributed uniformly in its envelope, 
and with limits on a dense central core amounting to > 0.02 by mass 
fraction (or about 6 earth masses maximum) (Guillot, Gauticr, and 
Hubbard, 1997). Models of Saturn, on the other hand, are distinctly 
nonsolar, with envelope enhancements of C, N, and O by about twice 
the Jovian factor, and a similar dense core of ^ 6 to 8 earth masses 
(Guillot et al. 1994). Uranus and Neptune models seem to resemble a 
larger version of the Saturn or Jupiter core with a thin (~ 0.1 by mass) 
envelope of hydrogen and helium (Hubbard, Podolak, and Stevenson 
1995). 

The traditional interpretation of this sequence, which is now open 

to revision in the light of recent detections of extrasolar giant planets 
and recent downward revisions of the Jupiter and Saturn core masses, 
is that the giant planets Jupiter and Saturn formed by the capture of 
nebular hydrogen and helium onto dense Uranus/Ncptune-like cores 
(Mizuno, Nakazawa, and Hayashi 1978), with subsequent significant 
accretion of icy cometesimals being responsible for enrichment of C, N, 
and 0-bcaring molecules in their envelopes. Presumably, the capture 
of nebular hydrogen onto proto-Uranus and proto-Neptune was less 
eflicient owing to lower nebular density at their orbital radii and slower 
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accretion of their dense icy cores. 

The detection of EGPs at very small orbital radii has strongly 
suggested the possibility of significant radial migration of giant planets 
during the nebular phase (Trilling et al. 1998). Recent observations 
(Terebey et al. 1998) may support the possibility of direct formation 
of giant planets without the necessity of an initial dense core to trigger 
hydrodynamic accretion of hydrogen-helium gas. 

B. Constraints on Cooling Mecheinisms 

The luminosity of a giant planet is determined by the heat radiated 
into space by means of depletion of its interior entropy as determined 
by the atmospheric boundary condition. The latter is strongly affected 
by metallicity. Moreover, the heat evolved can be strongly increased 
by a redistribution of hydrogen (with high specific entropy) toward 
shallower, lower-temperature regions of the planet, while denser com- 
ponents (with lower specific entropy) sink to deeper, hotter regions. 
Thus, atmospheric depletion of metals and helium may be accompa- 
nied by higher luminosity. It is this effect which is believed to be 
responsible for the anomalously high luminosity of Saturn (Stevenson 
& Salpeter 1977). 

In principle, the detection of EGPs, whose values of luminosity, age, 
mass, and atmospheric composition can be determined, will provide 
a test of this concept along with constraints on the hydrogen phase 
diagram. The latter predicts that helium separation may be important 
for objects with M < O.OOIM0, while possible segregation of other 
species via a first-order phase transition in hydrogen may occur in 
objects with masses M < O.OIMq. 
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FIGURE CAPTIONS 

Figure 1. Evolution of the luminosity (in L©) of solar-met allicity M dwarfs 
and substellar objects versus time (in years) after formation (from Burrows 
et al. 1997). The stars, "brown dwarfs" and "planets" are shown as solid, 
dashed, and dot-dashed curves, respectively. In this figure, we arbitrarily 
designate as "brown dwarfs" those objects that burn deuterium, while we 
designate those that do not as "planets." The masses in Mq label most 
of the curves, with the lowest three corresponding to the mass of Saturn, 
half the mass of Jupiter, and the mass of Jupiter. 

Figure 2. Approximate Mie scattering efficiency of various-sized enstatite 
dust grains. Submicron grains can produce large scattering and extinction 
opacity at optical wavelengths and yet remain essentially invisible at longer 
wavelengths. However, because they scatter almost conservatively in the 
optical, silicates are poor candidates for the grains responsible for the lower 
optical flux of G1229 B. 

Figure 3. Gravity dependence of the H2O features in the 1.9 to 2.1 micron 
region. The three synthetic spectra correspond to models which all have 
the bolomotric luminosity of G1229 B. The top curve shows the Tcff — 
870 K, \ogg = 4.5 spectrum, used as a reference spectrum in this figure. 
The other curves axe flux ratios of the Toh — 940 K, \ogg = 5 spectrum 
to the reference spectrum, and the ratio of the Tes = 1030 K, logg = 5.5 
spectrum to the reference spectrum, shown by the middle and bottom 
curves, respectively. All spectra have been shifted vertically for clarity. 
The synthetic spectra are shown at a spectral resolution of 1500. 

Figure 4. Model thermal emission spectra for model atmospheres {g — 
1000 msec^'^, TcH = 1000 K) of varying metallicity. These low-resolution 
spectra demonstrate that broad band near-IR flux measurements will not 
be able to distinguish differences in metallicity below about about 1/10 
solar. For larger metallicities, J band is most sensitive to increasing metal- 
licity. 

Figure 5. Geometric albedo spectra for Uranus, Jupiter, and a model 7- 
Jupiter mass planet with Tes = 400 K. The model includes no clouds, and 
is thus darker than either planet longward of about 0.6 /xm. Absorption by 

methane and water removes incident photons before they can be Rayleigh 
scattered in the model. In contrast, Jupiter's and Uranus' cloud decks 
reflect brightly in between strong methane absorption bands. Jupiter is 
much darker than either the model or Uranus at blue and UV wavelengths. 
Dark photochemical hazes are predominantly responsible for lowering the 
reflectivity of the planet in this spectral range. 
Figure 6. Model spectra for 12 Mj, T^s = 300 K planet orbiting at 1 A.U. 
from a G2V star. Solid line shows reflected flux received at 10 pc if planet's 
atmosphere is cloud-free. Dotted line demonstrates enhancement in flux 
when water clouds arc added to the model. Both models arc from Marley 
et al. (1998). Long dashed line gives thermal emission in the case of no 
clouds (Burrows et al. 1997). For this object, reflected ucar-IR flux begins 
to dominate thermal emission when clouds form as the object cools through 
about 400 K. 
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Evolution of luminosity with time for different masses 
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